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Abstract 

We derive an exact and explicit Kraus decomposition for the reduced density of a 
quantum system simultaneously interacting with time-dependent external fields and 
a chaotic environment of thermodynamic dimension. We test the accuracy of the 
Kraus decomposition against exact numerical results for a CNOT gate performed 
on two qubits of an (7V + 2)-qubit statically flawed isolated quantum computer. Here 
the N idle qubits comprise the finite environment. We obtain very good agreement 
even for small N. 
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1 Introduction 



A statically flawed isolated quantum computer (QC) [1,2] can exhibit a rich 
variety of dynamical behaviors [3, 4, 5], many of which are detrimental to its 
operation. Large unitary shifts have been observed which can destroy the 
fidelity over the span of a single gate[3,4,5]. Non-unitary effects such as internal 
decoherence and dissipation can also affect performance [4,5]. In some cases the 
QC can act similar to a kicked top [6] in the exponential decay regime [5], in 
which case options for error correction strategies [1,7] are quite limited. In 
other configurations, part of a QC can be employed to measure the strength 
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of residual two-body interactions [3], and this knowledge may prove useful for 
error correction. Idle qubits can also be manipulated to improve performance 
of the active qubits[4,5]. In fact such self-interacting quantum systems have 
an exceedingly complex dynamics we are only beginning to explore. 

The statically flawed isolated QC can be viewed as a subsystem-bath model 
where the subsystem is the active part of the QC (i.e., a set of qubits which 
are manipulated with classical external fields) and the bath consists of the idle 
qubits. The flaws consist of single body imperfections in qubits and residual 
two-body interactions. Similar subsystem-bath models arise in many proposed 
quantum technologies which employ condensed phase architectures. Quantum 
control of a chemical reaction in molecules [8, 9] in a hosting medium (e.g. 
gas phase, solution, surface, solid), requires coherent manipulation of a few 
degrees of freedom (the reaction coordinates) while they simultaneously in- 
teract with the rest of the degrees of freedom of the molecule(s) and the 
degrees of freedom of the environment. The degrees of freedom of the rest of 
the molecule(s) are usually strongly self-interacting, and hence are not well 
modeled by independent harmonic oscillators. These give rise to an intrinsic 
decoherence dynamics[10] like that in a QC. Similarly, the degrees of freedom 
of the hosting medium are generally strongly self-interacting and anharmonic. 
These give rise to external decoherence. Thus, coherent quantum control in 
general consists of a coherently manipulated subsystem interacting with a 
strongly self-interacting multi-component bath. 

At high temperatures there are a number of theories which can be employed 
to model such complex dynamics. Redfield theory [11,12] and its generaliza- 
tions [13,14] take into account internal bath dynamics via realistic molecular 
dynamics simulations. However, Redfield theory is notorious for its positivity 
violation, incorrect prediction of the long time limit, and for having divergent 
correlation functions for baths with discrete spectra. The theory of Bulgac et 
al[15] which represents the environment as an ensemble of random matrices 
is also useful for simulation of complex environmental dynamics. The semi- 
classical Wigner method [16] can also prove quite accurate for moderately high 
temperatures. 

At very low temperatures the available computational schemes are more lim- 
ited. An approximate master equation[4,17,18,19,20] which is a non-Markovian 
generalization of completely-positive-dynamical-semigroup theory [21], has been 
shown to be accurate for some subsystems interacting with chaotic baths [4, 20]. 
A Kraus decomposition [22] for independent subsystems interacting with chaotic 
thermodynamic baths has also been developed [23], and shown to be accurate 
numerically. Neither theory is exact, at least for finite baths, and so such meth- 
ods are best used in tandem. Where they are in agreement some confidence 
can presumably be placed in their predictions. 
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Here we extend this Kraus decomposition to subsystems which also interact 
with classical external fields. This is of course the most important case for 
coherent control problems. We also consider a more general class of subsystem- 
bath interactions than was addressed in [23]. Finally, to obtain the greatest 
generality we make no assumptions regarding the Hermiticity of the subsystem 
Hamiltonian, so that the method can be employed with complex absorbing 
potentials [24]. The extended chaotic Kraus decomposition (CKD) we derive 
automatically satisfies all required conservation laws, i.e., Hermiticity, norm 
conservation (for Hermitian subsystems only) and complete positivity for the 
reduced density. It is only exact in the limit of a bath of thermodynamic 
dimension. We test the accuracy of our CKD against exact numerical results 
of a realistic isolated QC model and we obtain very good agreement. 

Our test model [4,5] represents a two-qubit register performing a CNOT gate 
while interacting with neighboring idle qubits via static internal residual inter- 
actions. Statically flawed isolated QC models offer a prototypical example of 
a manipulated subsystem interacting with a self-interacting, possibly chaotic, 
environment wherein exact quantum dynamics can be readily obtained. These 
exact numerical results serve as a benchmark against which new theories of 
complex environments can be tested. Accordingly, we test the time-dependent 
extension of our CKD against these exact benchmark results [4,5] for a large 
number of configurations: eight different initial states, two types of error gen- 
erating interactions (bit-flip and phase errors), and a number of intra-bath 
couplings in the chaotic regime. We obtain very accurate results despite our 
relatively small bath dimension. 

Organization of this manuscript is as follows: We present the time dependent 
extension of the CKD in section 2. In section 3 we discuss our isolated QC 
model. In sections 4 and 5 we test the Kraus decomposition against exact 
numerical results. In section 6 we discuss our results. 



2 Extension of chaotic Kraus decomposition 



We will now derive a CKD which can be applied more widely than that devel- 
oped in [23]. The general structure of the argument is similar to that of [23], 
and where the two coincide we will refer the reader to the previous manuscript 
for details. We will begin by deriving an abstract general Kraus decomposition 
which makes no assumptions about the properties of the Hamiltonians beyond 
that the total Hamiltonian is of the form 

H(t)=H s (t) + Y,S l *B IM + H B (1) 
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where Hs(t) = H s + £(t) consists of the time-independent native system 
Hamiltonian Hs, and the time-dependent Hamiltonian representing external 
driving fields £(t). and are interaction Hamiltonians in system and 
bath degrees of freedom, and H B is a bath Hamiltonian. We will then use 
properties specific to chaotic thermodynamic baths to transform the Kraus 
decomposition into a more tractable sum. 

We consider a product initial condition 

p(0) =p s (0)®p B (0) (2) 



where ps(0) is an arbitrary initial subsystem state, and ps(0) is an initial bath 
state of canonical form, i.e. 

-Ej/k B T 

p B (o) = T,—Q—m\, (3) 



where Q = X^exp {—E k /k B T}. Here \j) denotes a complete bath eigenbasis, 
i.e. H B \k) = E k \k) and J2k \k)(k\ = 1b- Obviously this is not the most gen- 
eral initial condition [25], nor is it the most likely initial condition. However, 
the CKD requires such an uncorrelated initial state, and this is a distinct 
limitation of the method. However, in many studies an uncorrelated initial 
state is assumed for simplicity or in ignorance of the true initial state, and so 
the CKD may prove very useful in spite of this obvious limitation. Also, the 
CKD we derive can be applied to one type of initially correlated state, e.g. 
p(0) = |Vi(0))(V>i(0)| ®p B {$) + \M0))(M0)\®Pb(0), where p B (0) and p|(0) 
are different thermodynamic chaotic baths of potentially differing tempera- 
tures. In this case, the CKD is applied independently to each term, and the 
sum is then an exact representation of pit). Finally, it is also worth pointing 
out that other master equation based theories do not require such restrictions 
on the form of the initial states[17,18, 19,20]. 



2.1 General Kraus decomposition 



Given these initial conditions and Hamiltonians the time evolved density is 
p(t) = U(t)p(0)W(t) where the propagator U(t) is given by 

t 

U{t) = texp [-(i/H) j H(t')dt'} (4) 

o 

and T denotes the time ordering operator. 
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The exact reduced density of the subsystem at time t, f>s(t), can be expressed 
by as a partial trace over the bath degrees of freedom 

p s (t) = Tr B {U(t)p(0)U\t)}. (5) 



Performing the partial trace operation in the bath eigenbasis gives 

Ps(t) = J2(j\U(t) (ps(0) ® ^H^l) U\t)\j). (6) 



Defining the Kraus operators fcj,k(t) via 

t iJk {t) = y/pi(J\U{t)\k), (7) 

where p k = exp{—E k /kBT}/Q are the initial populations of the bath eigen- 
states, it follows that Eq. (6) can be written in the Kraus decomposition 
form [22] 

Ps(t) = Y,^At)psml k (t). (8) 

Note that if the subsystem Hamiltonian is Hermitian then it is always true 
that 

££;,*(*)4*(*) = 4 

3,k 

but Hamiltonians with absorbing potentials [24] do not satisfy this relation. 
This is no way limits the usefulness of the decomposition. 

2.2 Chaotic Kraus decomposition 



We next employ properties of the coupling matrix elements of a chaotic ther- 
modynamic bath to simplify the forms of the Kraus operators tCj t k(t) . First 
we rewrite the total Hamiltonian (1) in terms of the bath eigenvalues and 
eigenvectors, i.e., 

H = H s (t) + E ^ E B i k \j) (k\ + E Ej \j) (J I , (9) 

A 4 j,k j 

where B^ k = (j\B fJ \k) are the matrix elements of the bath coupling operator 
in the complete bath eigenbasis. 
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Open system dynamics is largely governed by the properties of the coupling 
operators and B^. If a bath is chaotic with N degrees of freedom, then it can 
be argued[23] that as a result of scaling the off-diagonal matrix elements of any 
operator B, in the bath eigenbasis, will vanish as N — > oo (i.e. (j\B\k) — > 
for j 7^ k). This result has been known for some time [26], but it has an 
important application here. Applied to our bath coupling operators it means 
that B^ k = for j '• ^ k in the limit of a thermodynamic bath. 

Using this property, the total Hamiltonian (9) simplifies to 

H = H s (t) + E(E W + E k)\k)(k\. (10) 

k |U 

Taylor expansion of (7) using the Magnus formula[27] then gives 



£j,k(t) = v^O'lfl + HA) J dh H(h) 

o 

t h 

+(-i/h) 2 f dhf dt 2 TH(h)H(t 2 ) 

o o 

t ti t 2 

+ {-i/hf J dhj dt 2 J dt 3 TH(h)H(t 2 )H(t 3 ) + ...]\k) (11) 



which, since H(t) is block diagonal in the bath eigenbasis, simplifies to 



z 

£j,k(t) = Vp~k[i + (-i/h) J dt, (j\H(h)\j) 

o 

t ti 

+(-i/h) 2 J dh J dt 2 fma^ma^j) 



t ti t 2 

\3 



H-i/W j dh J dt 2 J dt 3 T(j\H(h)\j)(j\H(t 2 )\j)(j\H(t 3 )\j) 



+ ..]5 hk . (12) 
Now, using the Magnus Taylor expansion [27] in reverse gives 



£j,k(t) = VPk~fexpl-^J dt'(H s (t') + £ SpB h * + E h ) \ 5 hk (13) 



6 



which in turn can be substituted into (8) to obtain the final form of the CKD. 
The double sum becomes a single sum and 

P5(*) = E&,*(*)&(0)Q >Jfc (f) (14) 



with iC kjk (t) given by (13). 

Eq. (14) is our final result. It extends our previous CKD [23] to time-dependent 
subsystem Hamiltonians, potentially non-Hermitian subsystem Hamiltonians, 
and more general subsystem-bath coupling operators. 



2.3 Numerical strategy 



While Eq. (14) is exact only for thermodynamic chaotic baths, our previous 
experience in [23] suggests that it may be accurate even for quite small baths. 
Thus, it is important to address the issue of how (14) can be employed in 
practice. 

First, note that the sum over k can probably be truncated at quite small 
values for very low temperature systems like QCs. Furthermore, if ps(0) = 
\ip(0)){ip(0)\ then one can define states 

\Mt))=te X p\-y dt'(H s (t')+j:S,B k /)\ \m) (15) 



o 



such that 



ps(t) = Y,Pk\Mt))(Mt)\ (16) 

k 



and 

d\Mt))/dt = -(t/h)[H s (t) + J2S,B k ^}\Mt)}. (17) 

Now, equations (17) can be solved exactly using standard Runge-Kutta[28] 
techniques. In the case where the initial state is not pure one could obviously 
find a similar set of Liouville-von Neumann equations, 

dp k s (t)/dt = -(i/h)[H s (t) + £ SpBf, p k s (t)l (18) 
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with p|(0) = ps(0), from which the reduced density can be constructed via 
Ps(t) = J2kPkf>s(t)- Once again these equations can be solved numerically 
with standard techniques. 



3 Test model 

Our test model represents a two-qubit register performing a CNOT gate while 
interacting with idle neighboring qubits via static residual interactions. We 
consider a two dimensional circuit and so the number of idle qubits iV is 
ten, which clearly does not constitute a bath of thermodynamic dimension. 
It is however chaotic when the magnitude of the residual interactions is suffi- 
ciently large. This can be verified using the nearest neighbor eigenvalue spacing 
distribution [29, 30] and using Loschmidt echo[31,32] calculations [23]. Note that 
we do not repeat these calculations here since the bath is the same in this and 
the previous study [23]. We will see that the CKD is accurate in the chaotic 
regime even though the bath dimension is small. 

We consider the following total Hamiltonian for the (iV + 2)-qubit isolated QC 
H(t)=H s (t) + H SB + H B . (19) 

Here Hs(t) represents the time-dependent control Hamiltonian that imple- 
ments the CNOT gate, H B is a bath Hamiltonian which represents the idle 
qubits, and Hsb is the interaction Hamiltonian which is responsible for the 
generation of errors. 

Elementary gate operations comprising the CNOT protocol are implemented 
using the control Hamiltonian [33] 

H S (t) = -\ E(fiT(*)*i + + j*(t)*&- (20) 

z 1=1 

In our CNOT implementation, we assume that the control Hamiltonian (20) is 
free of all imperfections, and that gates can be switched on and off by perfect 
square pulses. We also do not allow free system evolution. These conditions 
ensure a perfect implementation with maximum fidelity. Therefore, any er- 
ror that is observed in the CNOT implementation, regardless of its type (i.e. 
unitary or non-unitary), is due to the subsystem-bath interactions. The sub- 
system elementary gate Hamiltonians and their corresponding time intervals 
can be found in Table I. 
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Table 1 

The switching intervals and active Hamiltonians used to implement CNOT gate. 



Switching Intervals 

[70 = 0,71=^/(2^)] 

[t 1 ,t 2 = t 1 + it/(2B x )} 

[t 2 ,73 = 72 + 7T/(2£P)] 

[t 3 ,u = t 3 + V2tt/(2B z )} 
[T 4 ,r 5 = t 4 + it/(4J x )] 
[t 5 ,t 6 = t 5 + V2tt/(2B z )} 
[r 6 ,7v = r 6 + 7r/(2B*)] 
[t 7 ,t 8 = t 7 + it/(2B*)] 
[t$,t 9 = t$ + tt/{2B z )} 



Active Hamiltonian 

-\B*ol 

-\B*ol 
+\B*a 2 x 
+\B*al 



We consider two types of error generation, modeled by the interaction Hamil- 
tonian 

Hsb = (*i + d* a )± a , (21) 



where S a = I]^ 2 A^cr^ with index a G {x, z}. Henceforth, we sometimes refer 
to the errors generated by xx-type coupling as bit-flip errors and the errors 
generated by zz-type coupling as phase errors. 

While we assumed a perfect control Hamiltonian we did not make the same 
assumption for the bath Hamiltonian. We let the bath Hamiltonian include 
one and two-body static internal flaws. It has been argued in a number of 
studies [2] that one and two-body internal flaws are unavoidable features of 
multi-qubit QCs. Our bath Hamiltonian is then an iV-qubit generalization of 
the two-qubit control Hamiltonian (20) given by 



Hn = -- 



N+2 

E 

=3 



N+1 N+2 

(B*a x + B*&1) + £ £ X 

i=3 j=i+l 



x u x u x- 



(22) 



Here one-qubit parameters are sampled randomly and uniformly as Bf G 
[Bq — 5/2, Bq + 5/2] where Bq is the average value of the distribution, and 
5 is a detuning parameter. Two-qubit residual interactions and system-bath 
interactions are also sampled randomly and uniformly as J* J G [—J x , Jx] an d 
X l a G [—A, A], respectively. We will consider a number of J x values to explore 
the transition to strong chaos. 

We considered two sets of four initial system states ps(0) = IV'oXV'oI- The first 
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set consists of four standard basis states 

|Vo) e{|00),|01),|10>,|ll)}, (23) 

and the second set consists of four Bell states 

|Vo> G {(|00) ± |11»/ A (|01) ± |10»/V2}. (24) 

Thus, we have a set of benchmark tests which includes two types of subsystem- 
bath coupling, eight different initial states, and three J x values. In total, this 
set of tests is thus quite stringent. Reporting all of this data would however 
be problematic. Fortunately, there is very little state specificity in the data 
and so we can just report average quantities. That is, for each observable we 
average the results over the four standard initial states and the four Bell initial 
states, separately. 

3. 1 Quantitative error measures 

We quantified the extent of deviation from the ideal system evolution by two 
error measures: average purity and fidelity. Averages are taken over the four 
states of each set, i.e. over the four initial standard basis states, and over the 
four initial Bell states. 

The average purity is defined by 

= 7£Trs[pl(0] (25) 

The purity measures the degree of deviation from unitary dynamics and thus 
it quantifies non-unitary errors such as decoherence and dissipation. 

The average fidelity is defined by 

m= 1 jT,^s[ps(t)p l s deal (t)] (26) 

IY<o> 

where p l g eal (t) is the ideal subsystem evolution in the absence of interactions 
with the idle qubits. The fidelity measures how close the actual density stays 
to the ideal density in the course of the dynamics. Hence, it is sensitive to 
both unitary and non-unitary errors. Since the purity is insensitive to unitary 
effects, large deviations between purity and fidelity can be used as an indicator 
of unitary errors induced by the coherent shift process [3, 4, 5]. 
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The ideal value of the purity /fidelity is unity for pure initial conditions. 

We will also examine individual matrix elements, to see how the errors are 
manifested in the reduced density. 

3.2 Numerical parameters 

In our simulations we used the experimentally accessible control parameters 
for the charge-qubit QCs [33,34] for which B a = 1.00 e, J x = 0.05 e and 
kT = 0.25 e, in units of e = 200 mK. We assumed that the bath qubits 
only differ from the control qubits by static imperfections. Hence, we set the 
average value of the distribution to Bq = B a , and the detuning to S = 0.4 e. 
Our simulations included a number of two-qubit interactions. Here we only 
report our results for J x = 0.50, 1.00, 2.00, with units in e, for which we 
confirmed the onset of chaos in previous studies [3,4,23]. 

3.3 Exact numerical approach 

Exploiting the low temperature limit for the bath density, the time evolved 
subsystem reduced density can be expressed exactly as 

p s (t) = n jr Pn Tr B [\* n (t))(* n (t)\], (27) 

n=l 

where the populations of the bath are given by 

e -E n /k B T 

Pn = V n -9 p — Em/kgT ' ( 28 ) 

and where the states \^ n (t)) evolve from the Schrodinger equation with the 
total Hamiltonian (19) and initial conditions |\& n (0)) = \ipo) <8> \n). In all 
calculations n eig = 20 was sufficient for the given low temperature. We used 
a Lanczos algorithm[35] for exact diagonalization of the bath Hamiltonian 
and an eight-order variable stepsize Runge-Kutta code [28] for the numerical 
integrations. 

3.4 Chaotic Kraus decomposition approach 

Analytic solutions for the chaotic Kraus decomposition should be readily ob- 
tainable due to the low dimension of the CNOT system. However, we em- 
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ploy an alternative approach here by employing numerical solutions of the 
Schrddinger equation. This is more convenient due the number of initial states, 
couplings, and observables. Moreover, this approach is generally applicable for 
any subsystem. Since the initial states are all pure we used Eq. (17) and the 
sum (16). 

We used the n e i g = 20 low lying exact bath eigenstates |n) for J x = 0.50, 1.00, 2.00 
to calculate the B n,n required by the CKD. For xx and zz type coupling op- 
erators B n ' n = £™' n = (n|E a |n) where a G {x, z} stands for bit-flip and phase 
type couplings, respectively. 

The CKD requires exact bath eigenstates in the calculations of the B™ ,n terms. 
Exact diagonalization of the bath Hamiltonian should be easy to achieve 
for quite large bath dimensions by standard Lanczos matrix diagonalization 
routines[35], or more recent generalizations [36]. In cases where the bath di- 
mension is too large, e.g. a large coupled oscillator bath where exact diago- 
nalization is impossible, alternative approaches may be taken to calculate the 
_B"' n and E n . For example, hybrid quantum-semi-classical molecular dynamics 
simulations can prove very useful to calculate B™ ,n terms. Perhaps even the 
Wigner method[16] would suffice to calculate B^ n once the E n are known. 
These approaches will be investigated elsewhere. 



4 Results for average purity and fidelity 

We plot average purity V(t) for xx-type bit-flip coupling in Figure 1 and for 
zz-tjpe phase coupling in Figure 2 for three different values of intra-bath 
coupling, J x = 0.50, 1.00, 2.00. Decoherence and dissipation result in a pu- 
rity decay of less than 1 % over the course of the gate. The exact numerical 
results are represented by solid lines and Kraus results are represented by 
dotted lines, and each coupling value J x is assigned to the same color/line 
convention throughout. Results for standard basis states and Bell states are 
shown in subfigures (a) and (b), respectively. Switching times of elementary 
gate operations are also indicated by the grid lines. 

Figure 1 (a) shows excellent quantitative agreement between the exact and 
CKD for the most chaotic case of J x = 2. For J x — 1 the agreement is also 
quite good with errors in the purity of less than .05 %. For the least chaotic 
case of J x = .5 the discrepancy is on the order of .2 %. On the whole these 
results are surprisingly accurate for such a small bath. The results for Bell 
states in Figure 1 (b) are virtually indistinguishable from those of (a). 

Figure 2 (a) is quite different from Figure 1 (a). The overall exact decay of 
purity is comparable, but the CKD discrepancies for J x — 2 are on the order 



12 



(a) 

1.000 
0.998 
0.996 
0.994 
0.992 
0.990 
0.9S8 

1.000 
0.998 
0.996 
0.994 
0.992 
0.990 
0.988 
0.986 



Fig. 1. Exact numerical (solid lines) and Kraus (dotted lines) results for averaged 
purity V(t) in case of xx-type coupling for J x = 0.50 (black), J x = 1.00 (red), and 
J x = 2.00 (blue), (a) Standard basis states and (b) Bell states. 

of .1 %, while those for J x = 1 are about .2%. For J x = .5 the error is terrible. 
Figure 2 (b) for Bell states is again the same as (a) for standard states. 

We plot T{t) for xx-type coupling in Figure 3 and for zz-type coupling in 
Figure 4, for the same J x values. Here, we see perfect agreement between 
exact and CKD predictions for all states, all J x values, and all couplings. This 
is all the more remarkable when one notes that these errors are now very 
large. The fidelity in Figure 3 decays to 10 % of its initial value, while that 
in Figure 4 decays to 70 % of its initial value. The average fidelity of Figure 
3 for xx coupling shows no sensitivity to J x , while that in Figure 4 for zz 
coupling varies substantially with J x . The CKD captures both these effects. 
The large magnitude of the infidelity in both cases is due to a coherent shift 
of the subsystem[3,4,5]. Detailed discussions of the reasons for the different 
behaviors for the different coupling are given in Ref. [5]. 



5 Results for matrix elements of reduced density 

Purity and fidelity suffice for an overall identification of the magnitudes of 
non-unitary and unitary errors. The comparison of ideal and actual reduced 
density matrix elements, however, provides further valuable information on 
what actually goes wrong in an algorithm during open system dynamics. Here 
again it is impossible for us to present all of this data. As a generic repre- 
sentation of our results, we present a comparison of these matrix elements for 
two initial subsystem states; 1 11) is representative of the standard basis states, 
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Fig. 2. Exact numerical (solid lines) and Kraus (dotted lines) results for averaged 
purity V(t) in case of zz-type coupling for J x = 0.50 (black), J x = 1.00 (red), and 
J x = 2.00 (blue), (a) Standard basis states and (b) Bell states. 



(a) 




Fig. 3. Exact numerical (solid lines) and Kraus (dotted lines) results for averaged 
fidelity J-(t) in case of xx-type coupling for J x = 0.50 (black), J x = 1.00 (red), and 
J x = 2.00 (blue), (a) Standard basis states and (b) Bell states. 

and (| 00) + |ll})/\/2 is representative of Bell states. We consider both xx-type 
and zz- type couplings for J x = l.OOe. 

We compare the matrix elements for xx-type coupling in Figure 5 and zz-type 
coupling in Figure 6 for the initial state |11). In these and subsequent figures 
the error-free coherent time evolution is given by black solid lines. Exact open 
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Fig. 4. Exact numerical (solid lines) and Kraus (dotted lines) results for averaged 
fidelity T{t) in case of zz-type coupling for J x = 0.50 (black), J x = 1.00 (red), and 
J x = 2.00 (blue), (a) Standard basis states and (b) Bell states. 

dynamics time evolutions are denoted by solid green lines, and the CKD results 
are given by dashed red lines. 

Each subfigure in a figure represents a different matrix element. The specific 
matrix element plotted in each subfigure is as follows: in subfigure (a) we 
plot Poo(t) = (0|Tr 2 [/3s(t)]|0). Similarly, in subfigure (b) we plot Pn{t) = 
(l|Tr 2 [p s (t)]|l>, in subfigure (c) Re{p$(t)} = Re{(0|Tr 2 [p 5 (t)]|l)} and in 
subfigure (d) 1m{p$(t)} = Im{(0|Tr 2 [p s (t)]|l)}. Similarly, in subfi gure (e) we 
plot p$(t) = (0|Tn[^ s (t)]|0), in subfigure (f) pf}{t) = (1 |Tr x 1 1) , in sub- 
figure (g) Re{p$(t)} = Re{(0\Tri[pg(t)]\l)} , and in subfigure (h) lm{p$(t)} = 
ImiiO^lpami)}. 

The agreement between the exact and CKD results is excellent in all cases. 
The deviations from ideal free evolution are large in all cases, while the dis- 
crepancies between the exact and CKD results are basically negligible. The 
worst deviations are again seen in the zz coupling case but these are still very 
small. 

We show results for the matrix elements for the initial Bell state ( 1 00) + 
1 11)) /a/2 in Figure 7 for xx-type coupling and in Figure 8 for zz-type cou- 
pling. Here again the agreement between the exact and CKD predictions is 
astonishingly good. The only visible deviations occur for zz coupling. See Fig- 
ure 8 (a) and (c), for example, where there are some small deviations. 
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Fig. 5. Matrix elements of reduced density of first and second qubits for xx-type 
coupling. The initial state of system is |11) and the intra-bath coupling J x = 1.00 e. 




Fig. 6. Matrix elements of reduced density of first and second qubits for zz-type 
coupling. The initial state of system is |11) and the intra-bath coupling J x = 1.00 e. 



16 




Fig. 7. Matrix elements of reduced density of first and second qubits for xx-type 
coupling. The initial state of system is a Bell state of the form (|00) + |ll))/>/2 and 
the intra-bath coupling J x = 1.00 e. 




Fig. 8. Matrix elements of reduced density of first and second qubits for zz-type 
coupling. The initial state of system is a Bell state of the form (|00) + |ll»/\/2 and 
the intra-bath coupling J x = 1.00 e. 
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6 Summary and discussion 



We have considered a quantum system interacting simultaneously with time- 
dependent external driving fields and chaotic baths of thermodynamic dimen- 
sion. We derived a Kraus decomposition for the subsystem reduced density 
in an explicit and computationally tractable form which is exact for baths 
of thermodynamic dimension. We tested the accuracy of the Kraus decom- 
position against exact numerical results for a model system consisting of an 
isolated statically flawed QC performing a CNOT gate. We obtained quite 
accurate results for a large number of different configurations in spite of our 
small bath dimension. These promising results suggest that the Kraus de- 
composition can be a very useful and practical computational tool for low 
temperature simulations of open quantum systems. In particular, the CKD 
could prove quite useful for determining low temperature chemical reaction 
rates via the half-collision problem. 
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